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1. INTRODUCTION 

Asynchronous electrical drive system is found in enormous applications and is considered as well- 
established thank to the advancement in current loop control [1]. However, facing mechanical loads, speed 
control loop design takes a decisive role, especially when the non-rigidity and backlash of the motor-to-load 
connection is taken into account [2-5]. Flexible coupling phenomenon is studied in [6], the authors develop 
an elegant control scheme integrated with a state observer for estimating motor and load torques. Brock et al. 
[7] propose a filter accompanied with a neural adaptive controller to eliminate high resonance frequencies. 
The controlled system is analysed through extensive experiments proving the effectiveness of the proposed 
solution. For a class of multi-mass electrical drive systems where measurement is not always available, a 
system parameter identification method is presented in [8]. The authors employ pseudo random binary signal 
and chirp wave to excite necessary information. The approach might be difficult to apply in a wide range of 
practical systems. A 2DOF PI control is systematically designed in [9] for a two-mass drive coupled via a 
toothed belt. The tracking results are compared with 1DOF and 1DOF with feedforward controls confirm that 
the method exhibits robustness and load torque rejection. Due to flat property of the system, Thomsen et al. 
[10] formulate a flatness based control for a finite stiffness coupling drive system fed by an inverter with 
speed feedback signal. Backstepping control is developed for a multi-motor system which is vastly found in 
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rewinding systems in [11] suppressing flexible coupling phenomenon in linear speed and tension control 
system. In order to deal with nonlinear stiffness, the authors in [12] implement adaptive backstepping to 
reduce connecting shaft oscillations. However, system stiffness characteristic must be known to carry out the 
control, this requirement might be very challenging in practice. Recently, a active disturbance rejection 
control (ADRC) is considered as an alternative to classical PID control, an application of ADRC to two-mass 
system can be found in [13]. A fuzzy based control for a three-mass system is presented in [14, 15], where 
fuzzy term is selected to modify IPD control. Other control approaches can be found in [16-22]. 

In general, there is a lack of comprehensive study on how popular control methods react when 
applying to two-mass drive system operating in nominal range. In the paper, we carry out a hardware-in-the- 
loop based comparative analysis of speed responses using backstepping and flatness-based controls in 
nominal and field-weakening modes. Initially, we present the two-mass drive system in a view of ideal 
current loop. Subsequently, PID, backstepping, and flatness control design steps are given in details. Finally, 
hardware-in-the-loop results of three controls are compared and some important conclusions are drawn. 


2. TWO-MASS MODEL SYSTEM WITH IDEAL CURRENT LOOP RESPONSE 
Typical configuration of a two-mass system is shown in Figure 1. 


Load 
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Figure 1. Typical configuration of a two-mass system 


When operating in field weakening region, the model of the system can be given as [2]. 
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In which, isa; isq aredq components of the stator current; w,w, rotor and synchronous speeds, 
respectively; P,a W;,are rotor flux dq components;o is total leakage factor; T,is rotor time constant:uU.q, 
Usgare stator voltage dq components; L,is stator inductance, 1, zare motor and load speeds; #1, Pare 
motor and load angle accelerations; is rotor angle; dis the coupling shaft damping coefficient; cis shaft 
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stiffness. It can be seen that the system of dynamical (1) is of 7™ order. With the assumption of t perfectly 
designed current controller, the system model can be reduced as: 


dim 1. 1 
dt T T ™ 
S d F er a2 
So, a — —m 
= hh i A p A P2 A M (2) 
AP = 1 — Q2 
S d w d 1 
=Z = = ——m 
i i ee 


where: im = a The state (2) is of 4™ order, stator current isa is used to control the motor flux and isq is 
m 


dedicated to speed control. For control design purpose (2) is rewritten in the following state-space form: 
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3. SPEED CONTROL DESIGN 


3.1. PI control 
PI controller is designed according to the symmetric optimal standard, so we have the PI control 


structure shown as Figure 2. The design of PI for the system is well-established, readers can refer to [23, 24] 


for more details. 
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Figure 2. Block diagram of PI control algorithm of a two-mass system 
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3.2. Backstepping control 
For designing step, it is assumed that all feedback signals from the motor, the shaft and the load 


sides are available. As in Figure 2, the backstepping control is divided into two steps. In each step, a 
subsystem will be studied and controlled. Selecting a proper Lyapunov candidate function for each design 
step can guarantees the asymptotic stability of the corresponding subsystem. After completing all two steps, 
the asymptotic stability of the whole system can be achieved [25]. The speed control equation is given: 
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where My = ky(Wra/Lm)isqi Ko = (3/2)Zp(Lm/Lr) is motor torque. For ease of deploying 
control design, we define: 


X1 = Q2 

X> = Qı 

u = isq (6) 
y =X, = Q2 


where x, and xare state variables; u is the control input and y is the output (the load rotating speed). Since 
then the system is expressed in the strict feedback form as (7): 


X14 = (2) = fix x2) 


; då 
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In order to implement backstepping control we define: 


Zı = X1 — Xia = P2 — P2a 
Z2 = X2 — X2dka = P1 — Pidka (8) 


where X1g = Qzq1s the set value of the output variable and Xz2d4kķa = Ọidkats the virtual control. Selecting the 
Lyapunov candidate function as (9): 


1 
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Taking derivative of (9) results in: 
V; = Z4Z4 (10) 


Because of Øq = 0 so that (8) rendered as Z, = @>, substitution Z,from (5) into (10) yields: 
d. 1 d. 
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If V, = —k,z2 < 0, with Vk, > 0 then the stable condition is satisfied. Defining: 


c 1 d . 
By = 5 aA (13) 


Substitution of (13) into (12) results in: 
. J 
Pidka — os [k12Z, + Br] (14) 


The control specified in (14) ensures that the load speed tracks the desired value. Next, the control for the 
motor side is designed. It is noted thatx,g = @2zqis constant, so that from (5), (8), and (13) can be shortened as: 


. dq : d. 
4 =F = G2 = 7 Ot B (15) 
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Select the Lyapunov candidate function as (16): 

V> = V, + z? (16) 
Taking derivative of (16) gives: 

eo ie are (17) 
From (11) and (12), we can write as: 

V, = —kıż? + ACE — iaka) + Z222 (18) 
After some basic operations, it can be shown that: 

V, = —kıż? — kzz% + k3z% + AG — iaka) + 2222 (19) 
where we have defined @, — Ọidka = Z2, SO that: 

Vo =ekz =k Falon Ez + Z>) (20) 
At this point in order to renderV, < 0, the following conditions must be fullfied: 

k> > 0 and kZ, + ma +2 =0 (21) 
Conditions specified in (21) imply that: 

Ëa — Praca = —ke(G1 — Praka) — 55% 22) 
From (14), it indicates that: 

Pidka = -2 (kıż, + Br) and @ — Pidka = P1 + 2 (kız + Br) (23) 
Substituting (15) and (22) into (23) gives: 

Or = — (kr + ka)G1 — (Ë kika + 2a, — 2 (ky + ka) By — 2 By (24) 


To proceed to calculate the final expression of the control current, first of all, we substitution from the second 
equation of the system (5) into (24) we have: 
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From the expression (25) and using fundamental calculation, we will get the final expression of isq 1s: 
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Based on (26), the proposed control structure for the two-mass system driven by the induction motor is 
shown in Figure 3. 
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Figure 3. Block diagram of backstepping control algorithm of a two-mass system 


3.3. Flatness-based control (FBC) 
3.3.1. Feedforward control design 

Feedforward control can be considered as the inverse mathematical model of the control object and 
plays a decisive role in the FBC structure. From (5), the feedforward control can be derived as: 


di, 
une la =ln tT, F3 
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i sa TEA (27) 
Based on equation @, = — 2 Öz + a (my — m) [26], it can be deduced that: 
1 1 
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Assume that the load torque is not varying, then the load observer can be developed as: 
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dt 
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where My = Kylmisg. then from (28) and (29), it is straightforward to show that: 
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The error model is: 
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de 0 E, 0 
dea e T 
a = l; | A flm òm) (32) 
dt A 
The characteristic equation of (32) is: 
det[sI — A] = s? — ls — > = 0 (33) 
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By selecting l, andl, such that: 
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3.3.2. Design of the reference trajectories 
The equation x = P ( te “*) -reEeNju= Q( ay a) is also called the “inverse” 
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Trajectories of speed reference are given as: 
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3.3.3. Feedback control design 
To handle the flux and the speed of the motor, the well-known PI controller is chosen. The detail 
equations of the flux and speed controller are described by (37) and (38), respectively: 
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where: Ai, = isg — ie ; Aw, =w z — w3. Laplace transform of (38) provides the transfer function of 


error speed and flux: 


BONS) ps (39) 
Aisg(s) — s[JaJ287+(dJy+dJz)s +cJı+acJ2] 





Based on feedforward and feedback control design and trajectories of speed and flux reference the 
proposed control scheme for the two-mass system driven by the induction motor with perfect stator current 
controller is shown in Figure 4. 
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Figure 4. Block diagram flatness-based control of two-mass system 


4. RESULT AND DISCUSSION 

In this section, the proposed control strategy is verified by hardware-in-the-loop (HIL) simulation. 
The system hardware is simulated in real time on the HIL platform with a time step of lus that means very 
close to physical models, while the pulse width modulation carrier frequency is 5 kHz. Voltage and current 
controllers are implemented in DSP TMS320F2808. The simulation parameters are: rated power=2.2 kW; 
rated phase voltage=380/220 VRMS; rated frequency=50 Hz; d=0.015 Nm/rad; c=0.28 Nm/rad, Ji=7.455 
x10° Kgm?; J2=5.5913x10° Kgm?; with load torque is 1.5 Nm. 

Simulation scenario is: at t=Os, the magnetic current is formed. Then at t=0.5s, the motor starts to 
speed up to 1000 rpm. At t=2s, the reverse speed, the motor starts to reverse to -1000 rpm. In order to test the 
system ability, low speed range of | rpm is also applied in the simulation. This section provides evaluation of 
HIL results of speed response. The evaluation of difference speed control strategies for the two-mass system 
is shown in Figure 5 and Table 1. 


PI control Flatness-based control Backstepping method 


. 





+1000rpm 


-1000rpm 





(b) 
Figure 5. Speed response of Simulation/HIL platform at, (a) +1 rpm, (b) +1000 rpm 
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Table 1. HIL results of the speed response 


At the speed of +1 rpm At the speed of +1000 rpm 
Speed control Deadbeat PI Deadbeat Deadbeat Deadbeat PI Deadbeat Deadbeat 
Backstepping Flatness Backstepping Flatness 
Accelerating time (s) 50 50 0 50 50 0 
Speed settling time (s) 2.5 2.0 1.0 3.0 2.5 1.0 


It is found that the nonlinear control method based on the flatness-based control principle gives 
better performance than the PI and backstepping controls, as the referenced load speed match with motor 
speed after a short period of time, 0.3 seconds at (+) 1 rpm and 0.5 seconds at (+) 1000 rpm (at up and 
reverse speed), without overshoot. That implies the speed controller based on the flatness-based control can 
suppress resonance oscillation over the entire operating range of the IM motor. 

From the results, it can be seen that the two-mass system with flexible couplings when operating at 
the speed of (+) 1000 rpm and even at the low speed range (+) 1 rpm, the nonlinear control method, based on 
the flatness-based control, has proved its advantage by solving the problem of resonant vibration suppression 
at the spindle compared to the PI controller and backstepping. Besides of at times of transition with fast 
setting time, without overshoot and load speed match with motor speed. In the following simulation scenario, 
the motor is set in field-weakening region. The magnetizing current responses are shown in Figure 6 and 
Figure 7. The three control methods show stable operation not only in nominal but also field weakening 
range. However, it can be seen that the flatness-based control gives better transient response, 1.e., without 
overshoot, as seen in Figure 7, and shorter settling time (0.2s in comparison with 0.5s for PI control and 
backstepping method). 
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Figure 6. (a) Torque responses, (b) Speed responses 
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Figure 7. Controlled system performance, (a) Flux responses, (b) speed responses 
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Figure 7. Controlled system performance, (c) torque responses (continued) 


CONCLUSION 


In this paper, the two-mass system with flexible couplings comprising of an induction motor and a load 


is considered. At the motor side, the current response is assumed to be ideal leading to a simplified system 
model. PI control, Flatness-based control, and backstepping control are employed to solve flux and speed 
control problems of the system. The HIL based test shows better control performances of flatness control 
compared to the other two schemes. This conclusion implies possibility of practical uses of flatness control in 
dealing with systems with flexible couplings. In the future research, in order to give a concrete evaluation, 
robustness experiments against varying system parameter and load disturbances will be carried out. 
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